The Effect of Microwaves on Protein Structure: Molecular Dynamics Approach.

The impact of microwave (MW) irradiation on protein folding, potentially inciting misfolding, was investigated by employing molecular dynamics (MD) simulations. Twenty-nine proteins were subjected to MD simulations under equilibrium (300 K) and MW conditions, where the rotational temperature was elevated to 700 K. The utilized replacement model captures the microwave effects of δ- and γ-relaxation processes (frequency range of ∼300 MHz to ∼20 GHz). The results disclosed that MW heating incited a shift toward more compact protein conformations, as indicated by decreased root-mean-square deviations, root-mean-square fluctuations, head-to-tail distances, and radii of gyration. This compaction was attributed to the intensification of intramolecular electrostatic interactions and hydrogen bonds within the protein caused by MW-destabilized hydrogen bonds between the protein and solvent. The solvent-accessible surface area (SASA), particularly that of polar amino-acid residues, shrank under MW conditions, corresponding to a reduced polarity of the water solvent. However, MW irradiation produced no significant alterations in protein secondary structures; hence, MW heating was observed to primarily affect the protein tertiary structures.

Engineering A-type Dye-Decolorizing Peroxidases by Modification of a Conserved Glutamate Residue.

Dye-decolorizing peroxidases (DyPs) are recently identified microbial enzymes that have been used in several Biotechnology applications from wastewater treatment to lignin valorization. However, their properties and mechanism of action still have many open questions. Their heme-containing active site is buried by three conserved flexible loops with a putative role in modulating substrate access and enzyme catalysis. Here, we investigated the role of a conserved glutamate residue in stabilizing interactions in loop 2 of A-type DyPs. First, we did site saturation mutagenesis of this residue, replacing it with all possible amino acids in bacterial DyPs from Bacillus subtilis (BsDyP) and from Kitasatospora aureofaciens (KaDyP1), the latter being characterized here for the first time. We screened the resulting libraries of variants for activity towards ABTS and identified variants with increased catalytic efficiency. The selected variants were purified and characterized for activity and stability. We furthermore used Molecular Dynamics simulations to rationalize the increased catalytic efficiency and found that the main reason is the electron channeling becoming easier from surface-exposed tryptophans. Based on our findings, we also propose that this glutamate could work as a pH switch in the wild-type enzyme, preventing intracellular damage.

Molecular cloning, characterisation and molecular modelling of two novel T-synthases from mollusc origin.

The glycoprotein-N-acetylgalactosamine ß1,3-galactosyltransferase, known as T-synthase (EC 2.4.1.122), plays a crucial role in the synthesis of the T-antigen, which is the core 1 O-glycan structure. This enzyme transfers galactose from UDP-Gal to GalNAc-Ser/Thr. The T-antigen has significant functions in animal development, immune response, and recognition processes. Molluscs are a successful group of animals that inhabit various environments, such as freshwater, marine, and terrestrial habitats. They serve important roles in ecosystems as filter feeders and decomposers but can also be pests in agriculture and intermediate hosts for human and cattle parasites. The identification and characterization of novel carbohydrate active enzymes, such as T-synthase, can aid in the understanding of molluscan glycosylation abilities and their adaptation and survival abilities. Here, the T-synthase enzymes from the snail Pomacea canaliculata and the oyster Crassostrea gigas are identified, cloned, expressed, and characterized, with a focus on structural elucidation. The synthesized enzymes display core 1 ß1,3-galactosyltransferase activity using pNP-α-GalNAc as substrate and exhibit similar biochemical parameters as previously characterised T-synthases from other species. While the enzyme from C. gigas shares the same structural parameters with the other enzymes characterised so far, the T-synthase from P. canaliculata lacks the consensus sequence CCSD, which was previously considered indispensable.

Guidelines for Free-Energy Calculations Involving Charge Changes.

The Coulomb interactions in molecular simulations are inherently approximated due to the finite size of the molecular box sizes amenable to current-day compute power. Several methods exist for treating long-range electrostatic interactions, yet these approaches are subject to various finite-size-related artifacts. Lattice-sum methods are frequently used to approximate long-range interactions; however, these approaches also suffer from artifacts which become particularly pronounced for free-energy calculations that involve charge changes. The artifacts, however, also affect the sampling when plain simulations are performed, leading to a biased ensemble. Here, we investigate two previously described model systems to determine if artifacts continue to play a role when overall neutral boxes are considered, in the context of both free-energy calculations and sampling. We find that ensuring that no net-charge changes take place, while maintaining a neutral simulation box, may be sufficient provided that the simulation boxes are large enough. Addition of salt to the solution (when appropriate) can further alleviate the remaining artifacts in the sampling or the calculated free-energy differences. We provide practical guidelines to avoid finite-size artifacts.

Understanding the riddle of amine oxidase flavoenzyme reactivity on the stereoisomers of N-methyl-dopa and N-methyl-tyrosine.

Enzymes are usually stereospecific against chiral substrates, which is commonly accepted for the amine oxidase family of enzymes as well. However, the FsqB (fumisoquin biosynthesis gene B) enzyme that belongs to the family of sarcosine oxidase and oxidizes L-N-methyl-amino acids, shows surprising activity for both enantiomers of N-methyl-dopa. The aim of this study is to understand the mechanism behind this behavior. Primary docking experiments showed that tyrosine and aspartate residues (121 and 315 respectively) are located on the ceiling of the active site of FsqB and may play a role in fixing the N-methyl-dopa via its catechol moiety and allowing both stereoisomers of this substrate to be in close proximity of the N5 atom of the isoalloxazine ring of the cofactor. Three experimental approaches were used to prove this hypothesis which are: (1) studying the oxidative ability of the variants Y121F and D315A on N-methyl-dopa substrates in comparison with N-methyl-tyrosine substrates; (2) studying the FsqB WT and variants catalyzed biotransformation via high-performance liquid chromatography (HPLC); (3) molecular dynamics simulations to characterize the underlying mechanisms of the molecular recognition. First, we found that the chemical characteristics of the catechol moiety of N-methyl-dopa are important to explain the differences between N-methyl-dopa and N-methyl-tyrosine. Furthermore, we found that Y121 and D315 are specific in FsqB and not found in the model enzyme sarcosine oxidase. The on-bench and theoretical mutagenesis studies show that Y121 residue has a major role in fixing the N-methyl-dopa substrates close to the N5 atom of the isoalloxazine ring of the cofactor. Simultaneously, D315 has a supportive role in this mechanism. Jointly, the experimental and theoretical approaches help to solve the riddle of FsqB amine oxidase substrate specificity.

Size Matters: Free-Energy Calculations of Amino Acid Adsorption over Pristine Graphene.

There is still growing interest in graphene interactions with proteins, both for its possible biological applications and due to concerns over detrimental effects at the cellular level. As with any process involving proteins, an understanding of amino acid composition is desirable. In this work, we systematically studied the adsorption process of amino acids onto pristine graphene via rigorous free-energy calculations. We characterized the free energy, potential energy, and entropy of the adsorption of all proteinogenic amino acids. The energetic components were further separated into pair interaction contributions. A linear correlation was found between the free energy and the solvent accessible surface area change during adsorption (ΔSASAads) over pristine graphene and uncharged amino acids. Free energies over pristine graphene were compared with adsorption onto graphene oxide, finding an almost complete loss of the favorability of amino acid adsorption onto graphene. Finally, the correlation with ΔSASAads was used to successfully predict the free energy of adsorption of several penta-l-peptides in different structural states and sequences. Due to the relative ease of calculating the ΔSASAads compared to free-energy calculations, it could prove to be a cost-effective predictor of the free energy of adsorption for proteins onto nonpolar surfaces.

Amino Acid Residues Controlling Domain Interaction and Interdomain Electron Transfer in Cellobiose Dehydrogenase.

The function of cellobiose dehydrogenase (CDH) in biosensors, biofuel cells, and as a physiological redox partner of lytic polysaccharide monooxygenase (LPMO) is based on its role as an electron donor. Before donating electrons to LPMO or electrodes, an interdomain electron transfer from the catalytic FAD-containing dehydrogenase domain to the electron shuttling cytochrome domain of CDH is required. This study investigates the role of two crucial amino acids located at the dehydrogenase domain on domain interaction and interdomain electron transfer by structure-based engineering. The electron transfer kinetics of wild-type Myriococcum thermophilum CDH and its variants M309A, R698S, and M309A/R698S were analyzed by stopped-flow spectrophotometry and structural effects were studied by small-angle X-ray scattering. The data show that R698 is essential to pull the cytochrome domain close to the dehydrogenase domain and orient the heme propionate group towards the FAD, while M309 is an integral part of the electron transfer pathway - its mutation reducing the interdomain electron transfer 10-fold. Structural models and molecular dynamics simulations pinpoint the action of these two residues on the domain interaction and interdomain electron transfer.

Molecular modelling and site-directed mutagenesis provide insight into saccharide pyruvylation by the Paenibacillus alvei CsaB enzyme.

Pyruvylation is a biologically versatile but mechanistically unexplored saccharide modification. 4,6-Ketal pyruvylated N-acetylmannosamine within bacterial secondary cell wall polymers serves as a cell wall anchoring epitope for proteins possessing a terminal S-layer homology domain trimer. The pyruvyltransferase CsaB from Paenibacillus alvei served as a model to investigate the structural basis of the pyruvyltransfer reaction by a combination of molecular modelling and site-directed mutagenesis together with an enzyme assay using phosphoenolpyruvate (PEP; donor) and synthetic ß-D-ManNAc-(1 → 4)-α-D-GlcNAc-diphosphoryl-11-phenoxyundecyl (acceptor). CsaB protein structure modelling was done using Phyre2 and I-Tasser based on the partial crystal structure of the Schizosaccharomyces pombe pyruvyltransferase Pvg1p and by AlphaFold. The models informed the construction of twelve CsaB mutants targeted at plausible PEP and acceptor binding sites and KM and kcat values were determined to evaluate the mutants, indicating the importance of a loop region for catalysis. R148, H308 and K328 were found to be critical to PEP binding and insight into acceptor binding was obtained from an analysis of Y14 and F16 mutants, confirming the modelled binding sites and interactions predicted using Molecular Operating Environment. These data lay the basis for future mechanistic studies of saccharide pyruvylation as a novel target for interference with bacterial cell wall assembly.

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling.

Molecular dynamics simulations often struggle to obtain sufficient sampling to study complex molecular events due to high energy barriers separating the minima of interest. Multiple enhanced sampling techniques have been developed and improved over the years to tackle this issue. Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that works by adding a biasing potential, lifting the energy landscape up, and decreasing the height of its barriers. GaMD allows one to increase the sampling of events of interest without the need of a priori knowledge of the system or the relevant coordinates. All required acceleration parameters can be obtained from a previous search run. Upon its development, several improvements for the methodology have been proposed, among them selective GaMD in which the boosting potential is selectively applied to the region of interest. There are currently four selective GaMD methods that have shown promising results. However, all of these methods are constrained on the number, location, and scenarios in which this selective boosting potential can be applied to ligands, peptides, or protein-protein interactions. In this work, we showcase a GROMOS implementation of the GaMD methodology with a fully flexible selective GaMD approach that allows the user to define, in a straightforward way, multiple boosting potentials for as many regions as desired. We show and analyze the advantages of this flexible selective approach on two previously used test systems, the alanine dipeptide and the chignolin peptide, and extend these examples to study its applicability and potential to study conformational changes of glycans and glycosylated proteins.

Accelerated Enveloping Distribution Sampling to Probe the Presence of Water Molecules.

Determining the presence of water molecules at protein-ligand interfaces is still a challenging task in free-energy calculations. The inappropriate placement of water molecules results in the stabilization of wrong conformational orientations of the ligand. With classical alchemical perturbation methods, such as thermodynamic integration (TI), it is essential to know the amount of water molecules in the active site of the respective ligands. However, the resolution of the crystal structure and the correct assignment of the electron density do not always lead to a clear placement of water molecules. In this work, we apply the one-step perturbation method named accelerated enveloping distribution sampling (AEDS) to determine the presence of water molecules in the active site by probing them in a fast and straightforward way. Based on these results, we combined the AEDS method with standard TI to calculate accurate binding free energies in the presence of buried water molecules. The main idea is to perturb the water molecules with AEDS such that they are allowed to alternate between regular water molecules and non-interacting dummy particles while treating the ligand with TI over an alchemical pathway. We demonstrate the use of AEDS to probe the presence of water molecules for six different test systems. For one of these, previous calculations showed difficulties to reproduce the experimental binding free energies, and here, we use the combined TI-AEDS approach to tackle these issues.

Insights into the Binding Mode of Lipid A to the Anti-lipopolysaccharide Factor ALFPm3 from Penaeus monodon: An In Silico Study through MD Simulations.

The globally expanding threat of antibiotic resistance calls for the development of new strategies for abating Gram-negative bacterial infections. The use of extracorporeal blood cleansing devices with affinity sorbents to selectively capture bacterial lipopolysaccharide (LPS), which is the major constituent of Gram-negative bacterial outer membranes and the responsible agent for eliciting an exacerbated innate immune response in the host during infection, has received outstanding interest. For that purpose, molecules that bind tightly to LPS are required to functionalize the affinity sorbents. Particularly, anti-LPS factors (ALFs) are promising LPS-sequestrating molecules. Hence, in this work, molecular dynamics (MD) simulations are used to investigate the interaction mechanism and binding pose of the ALF isoform 3 from Penaeus monodon (ALFPm3), which is referred to as "AL3" for the sake of simplicity, and lipid A (LA, the component of LPS that represents its endotoxic principle). We concluded that hydrophobic interactions are responsible for AL3-LA binding and that LA binds to AL3 within the protein cavity, where it buries its aliphatic tails, whereas the negatively charged phosphate groups are exposed to the medium. AL3 residues that are key for its interaction with LA were identified, and their conservation in other ALFs (specifically Lys39 and Tyr49) was also analyzed. Additionally, based on the MD-derived results, we provide a picture of the possible AL3-LA interaction mechanism. Finally, an in vitro validation of the in silico predictions was performed. Overall, the insights gained from this work can guide the design of novel therapeutics for treating sepsis, since they may be significantly valuable for designing LPS-sequestrating molecules that could functionalize affinity sorbents to be used for extracorporeal blood detoxification.

Accelerated Enveloping Distribution Sampling (AEDS) Allows for Efficient Sampling of Orthogonal Degrees of Freedom.

One of the most challenging aspects in the molecular simulation of proteins is the study of slowly relaxing processes, such as loop rearrangements or ligands that adopt different conformations in the binding site. State-of-the-art methods used to calculate binding free energies rely on performing several short simulations (lambda steps), in which the ligand is slowly transformed into the endstates of interest. This makes capturing the slowly relaxing processes even more difficult, as they would need to be observed in all of the lambda steps. One attractive alternative is the use of a reference state capable of sampling all of the endstates of interest in a single simulation. However, the energy barriers between the states can be too high to overcome, thus hindering the sampling of all of the relevant conformations. Accelerated enveloping distribution sampling (AEDS) is a recently developed reference state technique that circumvents the high-energy-barrier challenge by adding a boosting potential that flattens the energy landscape without distorting the energy minima. In the present work, we apply AEDS to the well-studied benchmark system T4 lysozyme L99A. The T4 lysozyme L99A mutant contains a hydrophobic pocket in which there is a valine (valine 111), whose conformation influences the binding efficiencies of the different ligands. Incorrectly sampling the dihedral angle can lead to errors in predicted binding free energies of up to 16 kJ mol-1. This protein constitutes an ideal scenario to showcase the advantages and challenges when using AEDS in the presence of slow relaxing processes. We show that AEDS is able to successfully sample the relevant degrees of freedom, providing accurate binding free energies, without the need of previous information of the system in the form of collective variables. Additionally, we showcase the capabilities of AEDS to efficiently screen a set of ligands. These results represent a promising first step toward the development of free-energy methods that can respond to more intricate molecular events.

In silico investigation and surmounting of Lipopolysaccharide barrier in Gram-Negative Bacteria: How far has molecular dynamics Come?

Lipopolysaccharide (LPS), a main component of the outer membrane of Gram-negative bacteria, has crucial implications on both antibiotic resistance and the overstimulation of the host innate immune system. Fighting against these global concerns calls for the molecular understanding of the barrier function and immunostimulatory ability of LPS. Molecular dynamics (MD) simulations have become an invaluable tool for uncovering important findings in LPS research. While the reach of MD simulations for investigating the immunostimulatory ability of LPS has been already outlined, little attention has been paid to the role of MD simulations for exploring its barrier function and synthesis. Herein, we give an overview about the impact of MD simulations on gaining insight into the shield role and synthesis pathway of LPS, which have attracted considerable attention to discover molecules able to surmount antibiotic resistance, either circumventing LPS defenses or disrupting its synthesis. We specifically focus on the enhanced sampling and free energy calculation methods that have been combined with MD simulations to address such research. We also highlight the use of special-purpose MD supercomputers, the importance of appropriate LPS and ions parameterization to obtain reliable results, and the complementary views that MD and wet-lab experiments provide. Thereby, this work, which covers the last five years of research, apart from outlining the phenomena and strategies that are being explored, evidences the valuable insights that are gained by MD, which may be useful to advance antibiotic design, and what the prospects of this in silico method could be in LPS research.

A molecularly engineered, broad-spectrum anti-coronavirus lectin inhibits SARS-CoV-2 and MERS-CoV infection in vivo.

"Pan-coronavirus" antivirals targeting conserved viral components can be designed. Here, we show that the rationally engineered H84T-banana lectin (H84T-BanLec), which specifically recognizes high mannose found on viral proteins but seldom on healthy human cells, potently inhibits Middle East respiratory syndrome coronavirus (MERS-CoV), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (including Omicron), and other human-pathogenic coronaviruses at nanomolar concentrations. H84T-BanLec protects against MERS-CoV and SARS-CoV-2 infection in vivo. Importantly, intranasally and intraperitoneally administered H84T-BanLec are comparably effective. Mechanistic assays show that H84T-BanLec targets virus entry. High-speed atomic force microscopy depicts real-time multimolecular associations of H84T-BanLec dimers with the SARS-CoV-2 spike trimer. Single-molecule force spectroscopy demonstrates binding of H84T-BanLec to multiple SARS-CoV-2 spike mannose sites with high affinity and that H84T-BanLec competes with SARS-CoV-2 spike for binding to cellular ACE2. Modeling experiments identify distinct high-mannose glycans in spike recognized by H84T-BanLec. The multiple H84T-BanLec binding sites on spike likely account for the drug compound's broad-spectrum antiviral activity and the lack of resistant mutants.

Mutations Increasing Cofactor Affinity, Improve Stability and Activity of a Baeyer-Villiger Monooxygenase.

The typically low thermodynamic and kinetic stability of enzymes is a bottleneck for their application in industrial synthesis. Baeyer-Villiger monooxygenases, which oxidize ketones to lactones using aerial oxygen, among other activities, suffer particularly from these instabilities. Previous efforts in protein engineering have increased thermodynamic stability but at the price of decreased activity. Here, we solved this trade-off by introducing mutations in a cyclohexanone monooxygenase from Acinetobacter sp., guided by a combination of rational and structure-guided consensus approaches. We developed variants with improved activity (1.5- to 2.5-fold) and increased thermodynamic (+5 °C T m) and kinetic stability (8-fold). Our analysis revealed a crucial position in the cofactor binding domain, responsible for an 11-fold increase in affinity to the flavin cofactor, and explained using MD simulations. This gain in affinity was compatible with other mutations. While our study focused on a particular model enzyme, previous studies indicate that these findings are plausibly applicable to other BVMOs, and possibly to other flavin-dependent monooxygenases. These new design principles can inform the development of industrially robust, flavin-dependent biocatalysts for various oxidations.

Identification of Activating Mutations in the Transmembrane and Extracellular Domains of EGFR.

The epidermal growth factor receptor (EGFR) is frequently mutated in human cancer, most notably non-small-cell lung cancer and glioblastoma. While many frequently occurring EGFR mutations are known to confer constitutive EGFR activation, the situation is less clear for rarely detected variants. In fact, more than 1000 distinct EGFR mutations are listed in the Catalogue of Somatic Mutations in Cancer (COSMIC), but for most of them, the functional consequence is unknown. To identify additional, previously unknown activating mutations in EGFR, we screened a randomly mutated EGFR library for constitutive EGFR phosphorylation using a recently developed high-throughput approach termed PhosphoFlowSeq. Enrichment of the well-known activating mutations S768I, T790M, and L858R validated the experimental approach. Importantly, we also identified the activating mutations S442I and L658Q located in the extracellular and transmembrane domains of EGFR, respectively. To the best of our knowledge, neither S442I nor L658Q has been associated with an activating phenotype before. However, both have been detected in cancer samples. Interestingly, molecular dynamics (MD) simulations suggest that the L658Q mutation located in the hydrophobic transmembrane region forms intermolecular hydrogen bonds, thereby promoting EGFR dimerization and activation. Based on these findings, we screened the COSMIC database for additional hydrophilic mutations in the EGFR transmembrane region and indeed detected moderate constitutive activation of EGFR-G652R. Together, this study demonstrates that unbiased screening for activating mutations in EGFR not only yields well-established substitutions located in the kinase domain but also activating mutations in other regions of EGFR, including the extracellular and transmembrane domains.

CASPON platform technology: Ultrafast circularly permuted caspase-2 cleaves tagged fusion proteins before all 20 natural amino acids at the N-terminus.

Fusion protein technologies improve the expression and purification of recombinant proteins, but the removal of the tags involved requires specific proteases. The circularly permuted caspase-2 (cpCasp2) with its specific cleavage site, efficiently generates the untagged protein. While cleavage with cpCasp2 is possible before all 20 proteinogenic amino acids, cleavage before valine, leucine, isoleucine, aspartate and glutamate suffers from slow, and before proline extremely slow, turnover. To make the platform fusion protein process even more general such that any protein with an authentic N-terminus can be produced with high efficiency, the bacterial selection system PROFICS (PRotease Optimization via Fusion-Inhibited Carbamoyltransferase-based Selection) was used to evolve cpCasp2 into a variant with a catalytic turnover two orders of magnitude higher and the ability to cleave before any amino acid. The high specificity and the stability of the original circularly permuted protease was fully retained in this mutant, while the high manufacturability was mostly retained, albeit with decreased soluble titer. Four point-mutations are responsible for this change in activity, two of which are located in or near the binding pocket of the active site. This variant was named CASPON enzyme and is a major component of the CASPase-based fusiON (CASPON) platform technology. Applicability for the production of recombinant proteins was demonstrated by enzymatic removal of the CASPON tag from five model proteins. The CASPON tag enables high soluble expressions, affinity purification and good accessibility for cleavage. The five industry-relevant proteins of interest were FGF2, TNF, GH, GCSF and PTH.

Clinical grade ACE2 as a universal agent to block SARS-CoV-2 variants.

The recent emergence of multiple SARS-CoV-2 variants has caused considerable concern due to both reduced vaccine efficacy and escape from neutralizing antibody therapeutics. It is, therefore, paramount to develop therapeutic strategies that inhibit all known and future SARS-CoV-2 variants. Here, we report that all SARS-CoV-2 variants analyzed, including variants of concern (VOC) Alpha, Beta, Gamma, Delta, and Omicron, exhibit enhanced binding affinity to clinical grade and phase 2 tested recombinant human soluble ACE2 (APN01). Importantly, soluble ACE2 neutralized infection of VeroE6 cells and human lung epithelial cells by all current VOC strains with markedly enhanced potency when compared to reference SARS-CoV-2 isolates. Effective inhibition of infections with SARS-CoV-2 variants was validated and confirmed in two independent laboratories. These data show that SARS-CoV-2 variants that have emerged around the world, including current VOC and several variants of interest, can be inhibited by soluble ACE2, providing proof of principle of a pan-SARS-CoV-2 therapeutic.

Inhibition of the Peroxygenase Lytic Polysaccharide Monooxygenase by Carboxylic Acids and Amino Acids.

Lytic polysaccharide monooxygenases (LPMOs) are widely distributed in fungi, and catalyze the oxidative degradation of polysaccharides such as cellulose. Despite their name, LPMOs possess a dominant peroxygenase activity that is reflected in high turnover numbers but also causes deactivation. We report on the influence of small molecules and ions on the activity and stability of LPMO during catalysis. Turbidimetric and photometric assays were used to identify LPMO inhibitors and measure their inhibitory effect. Selected inhibitors were employed to study LPMO activity and stability during cellulose depolymerization by HPLC and turbidimetry. It was found that the fungal metabolic products oxalic acid and citric acid strongly reduce LPMO activity, but also protect the enzyme from deactivation. QM calculations showed that the copper atom in the catalytic site could be ligated by bi- or tridentate chelating compounds, which replace two water molecules. MD simulations and QM calculations show that the most likely inhibition pattern is the competition between the inhibitor and reducing agent in the oxidized Cu(II) state. A correlation between the complexation energy and the IC50 values demonstrates that small, bidentate molecules interact strongest with the catalytic site copper and could be used by the fungus as physiological effectors to regulate LPMO activity.

BuRNN: Buffer Region Neural Network Approach for Polarizable-Embedding Neural Network/Molecular Mechanics Simulations.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations have advanced the field of computational chemistry tremendously. However, they require the partitioning of a system into two different regions that are treated at different levels of theory, which can cause artifacts at the interface. Furthermore, they are still limited by high computational costs of quantum chemical calculations. In this work, we develop the buffer region neural network (BuRNN), an alternative approach to existing QM/MM schemes, which introduces a buffer region that experiences full electronic polarization by the inner QM region to minimize artifacts. The interactions between the QM and the buffer region are described by deep neural networks (NNs), which leads to the high computational efficiency of this hybrid NN/MM scheme while retaining quantum chemical accuracy. We demonstrate the BuRNN approach by performing NN/MM simulations of the hexa-aqua iron complex.